AD-UIS  931 


UtaCLASblFIEO 


F0RCI9N  TCGHNOtMY  OlV  MltMT-^ATTERSON  AFS  OM  '  F/O  21/a  ^ 

METHOD  OF  CALCULATINO  FOCKET  ENOINC  NOKlE  CONTOUF  FOF  OPTIMUM  — ETC(U) 
AU9  82  T  MINHUA 

F  TO- 1 D ( RS ) T-0087-82  m 


I 


FTD-ID(RS)T-0057-82 


EDITED  TRANSLATION 


FTD-ID(RS)T-0057-82  3  August  1982 

MICROFICHE  NR:  FTD-82-C-001054 

METHOD  OF  CALCULATING  ROCKET  ENGINE  NOZZLE 
CONTOUR  FOR  OPTIMUM  THRUST 

By:  Tian  Miohua 

English  pages:  5 

Source:  Gongcheng  Rewull  Xuebao  i Vol.  2, 

Nr.  3,  August  1981,  pp.  259-261 

Country  of  origin:  China 
Translated  by:  Randy  Dorsey 
Requester:  FTD/TQTA 

Approved  for  public  release;  distribution  unlimited. 


THIS  TRANSLATION  IS  A  Rf  NOtTION  OR  TNI  ORIOI* 
NAL  RORIION  TEXT  WITHOUT  ANT  ANALTTICAL  OR 
lOITORIAL  COMMENT.  STATEMENTS  OR  TNBORIES 
AOVOCATSOORIMRLIEOARE  THOSE  OR  TNE  SOURCE 
AND  DO  NOT  NECSSURILV  RERLRCT  TNI  POSITION 
OR  ORNHON  OR  THE  ROREION  TtCNNQLOOV  OL 
VISION. 


TRANSLATION  DIVISION 
ROREION  TECNNOLOeV  DIVISION 
WR.ARO.  OHM. 


D  -ID(R8)T-0057-62 


Date  ^  ^  19 


GRAPHICS  DISCLAIMER 


All  figures,  graphics,  tables,  equations,  etc.  merged 
into  this  translation  were  extracted  from  the  best 
quality  copy  available. 


j  Accession  For 
j  gra&i 

;  TAB 

'  Ui::-»jv::orinced 
!  dur,Mf  i.cat  ion^ 


pi  ::ti'ibut  ion/ 


Avnil-bili^y  Codes 
'  tad/or 


METHOD  OF  CALCULATING  ROCKET  ENGINE  NOZZLE  CONTOUR  FOR  OPTIMUM  THRUST* 
by  Tian  Minhua 

1.  Principles  of  Calculating  Nozzle  Contour  for  Optimum  Thrust 


This  work  considers  an  ideal  two-dimensional,  steady-state,  Irrot- 
atlonal,  Isentroplc  flow.  A  flow  field  inside  a  nozzle  is  axlsymniet- 
rlc  and  we  chose  to  use  a  cylindrical  coordinate  system.  In  Figure  1 
ABTE  represents  the  line  of  intersection  between  the  contour  surface 

of  the  nozzle  and  the  meridian  plane. 

The  converging  section  of  the  nozzle  AB 
is  a  circular  arc  wall  surface  with  a 
curvature  radius  of  R^.  The  initial 
expansion  section  BT  downstream  of  the 
throat  is  a  circular  arc  wall  surface 
with  a  curvature  radius  of  R2.  We  chose 
a  nozzle  contour  TE  which,  at  a  given 
nozzle  length  L  and  ambient  pressure  , 
will  cause  the  thrust  produced  by  the 
nozzle  to  reach  maximum.  We  called  this  contour  "the  nozzle  contour 
for  optimum  thrust”.  We  chose  a  point  C  on  the  nozzle  axis  and  consid¬ 
ered  a  control  surface  EC  passing  through  a  point  E  at  the  nozzle  out- 
let.  This  control  surface  EC  and  the  right  characteristic  curve  pass¬ 
ing  through  circular  arc  separation  point  T,  intersect  at  point  D  (see 
Fig.  1).  Engine  thrust  R  can  be  obtained  from  the  aerodynamic 


*  This  paper  was  given  in  April, 1980  at  the  Third  All-China  Scienti¬ 
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parameters  at  control  surface  EC. 


«  -  j' [  (A  -  ?.)  +  2.yrfy 

JC  L  S10  A  J 


and  ^  represent  flow  pressure,  density,  velocity,  angle 
of  flow,  and  the  angle  of  Incidence  between  the  control  surface  and 
the  *  axis,  respectively.) 

Since  point  C  Is  fixed  and  nozzle  length  Is  Invariant,  the  nozzle 
length  after  point  C  Is  also  Invariant,  1.  e.. 


Xg—Xc 


cig  tf>dy  « 


constant 


Moreover,  due  to  the  continuity  of  the  flow  rate,  the  flow  rate 
past  the  control  surface  must  be  equal  to  the  flow  rate  past  throat 
cross  section  6,  1.  e.. 


G  _  r  — 21  Zieydy  -  constant 

Jc  m 


The  problem  winds  up  as  that  of  finding  the  condition  below  which 
satisfies  expressions  (2)  and  (3)  and  causes  thrust  R  to  reach  peak 
value,  1.  e.,  causes  functional  expression  (1)  to  take  a  conditional 
peak  value.  Using  the  Lagrange  product  factor^'3,  the  conditional 
peak  value  problem  becomes  an  unconditional  peak  value.  Since  we  can' 
obtain  a  variational  method  which  causes  thrust  R  to  reach  maximum 
value,  the  flow  parameters  at  control  surface  EC  should  satisfy  each 
of  the  following  conditions:  (given  an  ambient  pressure  o ) 

(1)  At  outlet  end  point  E  we  have: 

dn2ffij  “•  (2/*iW4)ctgat  (4) 

(2)  Control  surface  ED  serves  as  the  left  characteristic  curve: 

dy/dx  —  tg(e  +  •)  . 

i 

d6  ___  ctgg  dv  rinffwiia  _  ^  , 
dx  X  dx  yeM(9  +  •) 


(3)  Th«  aerodynamic  parameters  at  ED  also  satisfy  the  two  foil* 
owing  relationships: 

M*eM(9  ■)/eaia  *  MacMCdc  —  •■)/eM«g  .  (6) 
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(l  +  5^  A/»)  '■  .in>9ig«  -  y,Wi(l  +  fl.iga,  ( 7  ) 

where 

M*-( — 1 — y 
V  *  -  1  +  2/MV 

(4)  Due  to  the  continuity  of  the  flow  rate  the  following  expres¬ 
sion  clearly  holds  true: 

!1  '**'  -  fil  ^  -  r  pr  -  *  ■  2^d,  ( 8 ) 

}o  ta{,d-^a)  Jp  eotC^  —  a^  '  ' 

(c^a  and  M  represent  the  ratio  of  specific  heats  of  the  gas,  the 
Mach  angle,  and  the  Mach  nxunber,  respectively.) 

2.  General  Description  of  the  Calculation  Process 

In  nozzle  design,  the  initial  conditions  which  are  usually  given 
are  nozzle  throat  radii  R]^  and  R2,  ratio  of  specific  heats  of  the  gas 
s,  and  noz.'i  e  expansion  area  ratio  in  order  to  determine  the  noz¬ 
zle  contour  for  optimum  thrust.  With  regard  to  rocket  engine  nozzles, 
in  order  to  boost  the  thrust  and  decrease  the  length  it  is  necessary 
to  employ  a  relatively  small  throat  wall  curvature  radius.  We  employ¬ 
ed  the  conforihal  curve  coordinate  method  in  m  to  calculate  the  para¬ 
meters  at  Initial  transonic  line  GK  of  the  nozzle  throat  (Fig.  1) . 

This  method  of  calculation  is  simple,  its  accuracy  is  fairly  good, 
and  it  is  suitable  for  cases  where  the  throat  wall  curvature  radius 
is  relatively  small.  The  entire  supersonic  flow  field  was  calculated 
and  we  msde  use  of  a  characteristic  curve  graph As  to  the  left 
and  right  characteristic  curve  equations,  we  employed  difference  equa¬ 
tions  with  two-dimensional  accuracy,  1.  e., 

(The  parameters  with  represent  mean  values  of  a*) . 

The  following  relates  the  method  for  finding  the  nossle  contour 
for  optimum  thrust  TB  (see  Pig.  1) .  The  flow  parameters  at  right 
characteristic  curve  T'H',  which  runs  from  the  wall  to  the  axis,  were 
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determined  one  by  one  according  to  the  characteristic  curve  method. 

We  first*  assume  JWj’*  and  can  obtain  by  expression  (4).  Then  we 
look  for  a  point  D  on  T'H'  which  will  allow  expression  (6)  to  be  sat¬ 
isfied.  If  we  are  unable  to  find  point  D  on  T'H’,  then  we  can  again 
find  the  next  characteristic  curve.  When  we  find  point  D,  we  can  de¬ 
termine  by  expression  (7)  and  then  can  determine  the  pareuneters  of 
the  various  points  on  ED  by  expressions  (5-7)  and  at  the  seune  time 
use  expression  (8)  to  decide  whether  or  not  T*  Is  separation  point  T. 
If  expression  (8)  holds  true,  then  we  find  the  next  right  character¬ 
istic  curve  T'H'  and  point  D,  until  expression  (8)  Is  satisfied.  This 
time  T'H'  will  be  the  TH  right  characteristic  curve  which  is  called 
for.  Then  At*  is  determined  on  the  basis  of  y]^>-and  compared  with  the 
given  A,.  Normally  .4?*  will  not  be  equal  to  A,.  Using  the  "shooting" 
method  we  can  select  an  appropriate  WJ,*’,  repeatedly  using  the  above 
described  method  to  obtain  jfi*' until  it  is  sufficiently  close  to  A,. 

Then  once  more  using  the  parameters  at  right  characteristic  curve 
TO  and  control  surface  ED  as  the  initial  conditions,  on  the  basis  of 
wall  TE  serving  as  the  line  of  flow,  the  flow  field  parameters  of  tri¬ 
angular  region  TDE  (Fig.  1)  can  be  calculated  using  the  characteristic 
curve  method,  thereby  obtaining  the  nozzle  contour  for  optimum  thrust. 
TE. 

3.  Calculation  Results  and  Discussion 

Table  1  below  lists  the  separation  point  (point  T)  parameters  and 
thA  outlet  end  point  (point  E)  parameters  for  the  nozzle  wall  contour 
for  optimum  thrust  at  different  expansion  area  ratios  Ag  and  throat 
divergent  section  radii  R2.  In  the  table,  L  is  nozzle  length  and  c* 

Is  the  nozzle  thrust  coefficient. 

By  the  calculated  results  It  can  be  seen  that  when  the  area  ratio 
Is  fixed,  the  smaller  the  nozzle  throat  divergent  section  radius  R2» 
the  smaller  the  nozzle  length  L.  In  order  to  reduce  nozzle  length  we 
can  appropriately  reduce  R2*  But  R2  cannot  be  excessively  small  since 
the  smaller  R2  the  larger  the  expansion  angle  $,  at  the  initial  sep¬ 
aration  point.  When  R2  ia  too  small  the  throat  flow  can  expand  too 
quickly .causing  flow  losses  to  increase.  Consequently,  it  can  be 
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determined  by  tests  down  to  what  value  R2  will  be  more  suitable, 
der  normal  conditions,  the  throat  wall  curvature  radius  can  take 
Jti  <  2  and  0.5  <  X,  <  1. 
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Combining  the  eharncteriptic  method  with  variation  principle  and  applying  it  to  su* 
peraonic  flow  fieldp  a  eale.ulating  method  is  presented  for  the  optimum-thrust-nottle 
contour  design.  This  method  can  be  used  to  design  nossle  contours  with  various 
throat  vurvature  radii  and  expantion  area  ratios.  It  is  especially  suitable  f<ir  the 
d<*sifm  of  nosfJe  contour  with  large  expantlon  area  ratio. 


